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_L ■ Abstract 

a 

We study sums of directed paths on a hierarchical lattice where each bond has 



either a positive or negative sign with a probability p. Such path sums J have 
been used to model interference effects by hopping electrons in the strongly 



pg \ localized regime. The advantage of hierarchical lattices is that they include 

o, 

0^ \ path crossings, ignored by mean field approaches, while still permitting an- 

alytical treatment. Here, we perform a scaling analysis of the controversial 

a ■ 

"sign transition" using Monte Carlo sampling, and conclude that the transi- 
a , tion exists and is second order. Furthermore, we make use of exact moment 

o 

recursion relations to find that the moments ( J n ) always determine, uniquely, 
the probability distribution P{J). We also derive, exactly, the moment behav- 
ior as a function of p in the thermodynamic limit. Extrapolations (n — ► 0) to 
obtain (In J) for odd and even moments yield a new signal for the transition 
that coincides with Monte Carlo simulations. Analysis of high moments yield 
interesting "solitonic" structures that propagate as a function of p. Finally, 
we derive the exact probability distribution for path sums J up to length 
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L = 64 for all sign probabilities. 
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I. INTRODUCTION 

Sums of directed paths are present in numerous models of disordered systems. Polymer 
configurations in a disordered matrix, dynamics of interfaces grown by deposition Jl|] and 
Feynman path sums for electron hopping between impurities [^],[| are only a few of the 
relevant examples. In this paper, we focus on the latter example involving a model first 
introduced by Nguyen, Spivak and Shklovskii (NSS), for interference effects in the strongly 
localized regime M. 

In the directed path sign model one studies the sum of all possible directed paths between 
two sites on a lattice. On each lattice bond, one places a random sign with probability 
p. Each directed path evolved is then computed by multiplying the values of the bonds it 
crosses. Finally the sum J of all paths is obtained. The proponents of the model || obtained, 
numerically for small systems, that a second order transition occurred at p c ~ 0.05 between 
a phase with preferential sign (for the path sum J), and a phase with no preferential sign. 
NSS also offered appealing arguments based on the behavior of 5 J/ (J). Presumably, such 
a parameter grows exponentially above the transition while it goes to zero below p c . The 
physical relevance of this transition lies in the fact that it may signal the change between 
Aharonov-Bohm oscillations of period hc/e and those of hc/2e || in the context of hopping 
conduction. 

The NSS argument was later contended by Shapir and Wang f| arguing that correlations 
between paths implied that SJ/(J) does not necessarily go to zero for any p. Subsequently, 
Wang et al [f| used an exact enumeration scheme to probe the transition for small lattices 
of maximum size L = 9. The work found no evidence of a transition above negative sign 
probability p = 0.02. Such conclusions were supported by Zhao et al on the basis of 
numerics, for large square lattices, where it was assured that the transition did not exist 
above p = 0.025 in two dimensions. Nevertheless, the decay of the order parameter AP 
as a function of system size was found to be anomalously slow for finite p (see also ||). 
Therefrom, more recently Spivak, Feng and Zeng M discussed numerical results that suggest 



a finite jump in the order parameter indicating a first order transition for the sign problem. 
The authors also imply that the moments ( J n ) increase faster than n\ as n — ► oo indicating 
there is no unique relation between (J n ) and the probability distribution P(J). This is an 
important point since the moments, in such a case, may not contain information about the 
transition. Finally, in a recent paper by Nguyen and Gamietea || , a renewed extensive study 
of the parameter 5J/(J) proves that, at least according to such a parameter, no transition 
exists; only a strong crossover from logarithmic to exponential behavior is observed. 

Besides the numerical approaches, mean field type approximations by Obukhov || point 



to a second order transition for dimensions d > 4. Furthermore, Derrida and Cook |T0[ also 
took up the problem, analytically, using a sparse matrix approach. They generalized the 
model to random phases, which includes random signs as a special case. Their approach is 
mean field in nature and results in a phase diagram where the sign transition is of second 
order |IT| (see also fl2||). Nevertheless, mean field results may not apply to low dimension 



due to the importance of path crossings ||. 

Here we address the following issues: i) What is the order of the sign transition through 
scaling analysis of the order parameter proposed, ii) do moments of the path sums determine 
the probability distribution uniquely? and iii) what is the exact behavior of the parameter 
SJ/{J) above and below the transition?. A new perspective will be gained by using hierar- 
chical lattice: Such lattices, while still amenable to analytical manipulation, include crucial 
path correlation effects absent in mean field. 

The paper is organized as follows: Section II discusses the sign model and describes 
hierarchical lattices. In section III we perform detailed Monte Carlo simulations, close to 
the transition, for systems of up to size L = 512. A scaling analysis is performed for the 
order parameter AP = P(J > 0) — P(J < 0) to distinguish between first and second order 
transitions. In section IV we study the moments (J n ) exactly, using moment recursion 



relations [|H|. We find that moments determine the distribution uniquely according to 
Carlemans theorem, and find possible indications of a phase transition from odd and even 
moment extrapolations to n = 0. In this section we also discuss the high moment behavior, 



unveiling interesting structures as a function of the sign probability p. Subsequently, we 
probe the parameter 5 J / (J) exactly showing, for the first time, its unambiguous crossover 
between exponential and logarithmic behavior. In section V we obtain the exact probability 
distribution for lattice sizes L = 16 and sample the distribution for up to L = 64 as a 
function of p. We end with the conclusions and a discussion of the mapping of the moments 
to an n-body partition function in one dimension as a continuum model that might aid in 
explaining the curious high moment behavior. 

II. THE SIGN MODEL 

Imagine two reference points on a lattice between which one would like to evolve all 
possible directed paths and compute a "partition function" 

</ = Er, (i) 

i 

where Tj represents each individual contributing path. By directed it is meant that paths 
always propagate in the forward direction without loops or overhangs. The random medium 



in which these paths evolve can be represented by assigning local weights ]14[ on the bonds 
or sites that are picked up by the paths as they wander to their final destination. Such 
a model has been used as a paradigm simulating, for example, a coarse-grained polymer 
or interface wandering in a random matrix with locally favorable energy minima ||14|| . The 
model is interesting because it yields anomalous lateral wandering and energy exponents for 
the interface/polymer as compared to those generated by simple diffusion, signaling a new 
disorder induced universality class in (1 + 1) dimensions. 

Another application, in an entirely different field, is in the context of Variable Range 
Hopping [[HJ, a mechanism for conduction in insulators. In this context, one also needs 
to sum over Feynman paths to compute the transition probability, between impurities, of 
current bearing electrons. The Feynman paths, in this case, are directed because they 
are tunneling paths. Any elongation of the latter, in the form of loops or overhangs, is 



exponentially less probable. For further justification of the model we refer the reader to 
the review in reference ||. NSS studied such tunneling processes proposing a directed path 
model where the local weights are random signs 0. In such a model, the path Tj is a product 
of the signs it picks up en route to the final site. Writing eq.(jT|) more explicitly 

J = Y,IlVi, (2) 

r\ i 

where rji is a random sign according to the distribution P{rj) = p8(j] — 1) + (1 — p)S(r) + 1). 
The probability p in the NSS model emulates the relative abundance of levels above and 
below the Fermi energy . This model has been very successful in explaining qualitative and 
quantitative features of conduction in the strongly localized regime. In particular, intriguing 
interference effects producing a characteristic periodicity of magnetic field oscillations [16| 



and changes in the localization length due to non-local effects 0,0. In spite of the seemingly 



different nature of disorder in the NSS model, replica arguments and numerics have shown 
that it belongs to the same universality class of directed polymers with positive weights 
PJT8[], at least for p close to 1/2. 



We have taken up the sign model on hierarchical lattices as mentioned in the intro- 
duction. A hierarchical lattice is a recursive structure built by repeating a chosen motif 
T9fl . Depending on the latter motif, one can build integer dimensional objects emulating an 



Euclidean lattice. For this work we chose the Berker lattice or diamond. Such motif (see 
Fig. |l|) has the parameter b corresponding to the number of branches between the initial 
% and final / points. The lattice size is related to the recursion order m as L = 2 m_1 i.e. 
the number of bonds on any directed path between i and /. The number of bonds on the 
lattice (or mass) is given by M = (2b) m ~ 1 , so that the effective dimension of the lattice 
is d e ff = 1 + |^| . In this work we will use 6 = 2 except if otherwise stated. Qualitative 
features of critical behavior of many statistical models are correctly reproduced on such 
structures with no unphysical effects. In fact, mapping to hierarchical lattices is the basis of 
the Migdal-Kadanoff renormalization procedure, of frequent use in critical phenomena. As 
noted before, an important feature of hierarchical lattices over either Bethe lattices/mean 



field approaches is that path intersections are taken into account. Thus, we expect that 
the resulting simulations will be more faithful to low dimensional behavior. In fact, we 
will present, in section IV, further evidence of the adequacy of hierarchical lattices making 
contact with known recent results on the sign transition. 

III. SIGN PHASE TRANSITION IN TWO DIMENSIONS 

In this section, we have undertaken Monte Carlo simulations on hierarchical lattices to 
check for scaling properties. Paradoxically, scaling has only been discussed once before in 
connection with the transition |L8| and it is a primary tool to assess its nature. It will be 
especially useful to clearly distinguish between first order and second order transitions. 

Hierarchical lattices were generated to L = 512 or order ten. Averages were taken over 
twenty thousand realizations of disorder for a series of p values between and 0.5. As the 
size of the system increases more detailed data was collected close to the transition regime 
0.05 < p < 0.1. Figure |2] shows Monte Carlo data for the order parameter AP as a function 
of p. A definite plateau at AP = 1 develops as L increases for low p, signaling a definite 
change in the order parameter (positively signed paths dominate). 

For the proposed order parameter we should expect the scaling form AP = f{{p — 
Pc)L l / u ). Figure ^| shows a good collapse for the same data as the previous figure. As the 
order parameter is always between zero and one, we only need to find p c and the correlation 
length exponent v. For the hypothetical transition we find the values p c = 0.071 ±0.001 and 
v = 1.85 ± 0.07 (1/V = 0.54). The latter exponent is very different from that of percolation 
on these lattices [y = ln2/ (ln2 + ln(3 — V5)) = 1.63529..; so the role of percolation if 
any, is not apparent. If the transition were first order the exponent 1/V would be the 



dimensionality of the system d |20| . The non trivial scaling found can also be seen by taking 
the derivative of the order parameter and plotting its maximum as a function of the system 
size. These criteria rule out a first order transition. 

We have also monitored the evolution of p c (L) with size. The specific value of p c {L) 



was found from the peak values of the derivative of the order parameter AP. The resulting 
values are plotted in Fig. |] where, within error bars, the values of 1/v and p c (oo) are 
confirmed. Summarizing, scaling is very good around p c = 0.07 and does not correspond 
to the scaling of a first order transition. Furthermore, there is no sign of a discontinuity 
in the order parameter as suggested in ref. JF]. We thus conclude that, on hierarchical 
lattices, the transition exists and is second order as mean field predicts. These conclusions 



are in agreement with work by Roux and Coniglio |L8| on hierarchical lattices. There, they 
analyzed the variable cti = (nf — n~) where n i is the fraction of positive (negative) paths 
arriving at site i, and they suggest a clear positive a phase. The order of the transition for 
hierarchical lattices is not analyzed in detail in their paper. Nevertheless, they noted undue 
emphasis of hierarchical lattices on the a = result, and the possible impact of this on the 
scaling properties of various quantities. We will come back to such observations, briefly, in 
section V. 

IV. MOMENT RECURSION RELATIONS 

A statistic we can probe exactly on hierarchical lattices are the moments of the probabil- 
ity distribution. This is possible because of recursion relations derived by Cook and Derrida 
|2lfl and generalized to arbitrary moment and hierarchical order (system size) by Medina 



and Kardar |13|] . The recursion relation for b = 2 is the following 



n\ 



<-Cfi> = E ,„ \,„, [WH(-C~ S )] 2 , (3) 



__ {n-sy.s\ 



where n is the moment number and m is the hierarchical lattice order. This expression 
is readily generalized to other integer b by changing the binomial factor to a multinomial 
and including the additional branches. Hence, one can emulate higher dimensional networks. 
The simple form of this recursion permits, given the local moments at order one, to compute 
moments to any given lattice size. Appropriate programming of the recursion relations, with 
arbitrary precision computations, is linear in time with lattice order. 
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The behavior of the moments for the sign model is extremely rich as we shall see in 



the following. As found in ref. [[U|, after a few hierarchical orders, the values ln({J^ l ))/L 
converge rapidly to a limiting form as a function of n. Such limiting form is important 
because it also signals the convergence to a unique limiting distribution, at least, if moments 



do not grow faster than n\ |)22||. The asymptotic form of the moments can be obtained for 



p = UJ21 



*£>-(l-i)«. W 

that is, moments grow exponentially with n for p = 0. Nevertheless, for < p < 1/2, 
lower moments grow slightly faster than exponential (exp(n Q ) with 1 < a < 2), gradually 
converging to exponential growth for larger moments. The latter implies, according to the 
condition 

oo 

j2(J 2n r i/2n = oo, (5) 

n=0 

that the moments determine the distribution uniquely. There are various forms of such a 
theorem, but the above is the strongest version due to Carleman |2"2fl. If one substitutes 



above, (J 2n ) ~ exp(2n), -our asymptotic result- the criterion is satisfied. Even if (J 2n ) 
grows slightly faster i.e. exp(2nln2n) ~ (2n)\ the above sum diverges because J2 n V n = °°- 
Any faster growth would violate eq.(|5]), factorial growth being the borderline case. That the 
moments (J n ) satisfy eq.(|5]) is one of our central results. In Fig. [| we show a sequence of 
moments as a function of the moment number n. The different curves, starting from below, 
represent hierarchical orders one thru nine (sizes L = 2 thru L = 256). One readily notes 
convergence to a definite law. The inset shows a comparison between the growth of n\ and 
that of moments for the particular case of p = 0.1. The asymptotic behavior is already 
reached at L = 128, larger sizes falling on the same curve. 

For values close to p = 1/2, the moment sequence has a characteristic sawtooth shape, 
where even moments are at the crests and the odd at the troughs. Such structure is not 
a finite size effect. We have checked this for up to L = 2 20 on the hierarchical lattice. As 



p — ► 1/2 all the odd moments go to zero while the even remain finite as expected. On the 
other hand, as p is reduced the sawtooth disappears, first for the higher moments and then 
for the lower. In this respect there appears to be a phase transition for each moment at 
different values of p, in a way reminiscent of that discussed by Cook and Derrida |21| (in 
their case as a function of 'temperature'). The transition for the first two moments occurs 
close to p = 0.075 which is close to that found from Monte Carlo simulations in the previous 
section. On this basis it is plausible that the disappearance of sawtooth shape is related to 
the transition. 

Figure |6| shows a set of curves for d ln(( J^) / L) / dn, and various values of p, L = 2 17 and 
up to n = 100. The last six orders of the hierarchical lattice collapse on the same curve 
indicating we have achieved asymptotics. For the highest value of p one notes the sawtooth 
behavior, while it disappears for all moments below p = 0.1. Nevertheless, additional 
structure is observed at moments beyond n = 40 for p = 0.075 and p — 0.1, where a shoulder 
develops and moves towards larger n values as p increases, undeformed, in a solitonic manner. 
Although the analysis of these structures is out of the scope of this paper, it is interesting to 
analyze it in the light of a mapping to a one dimensional many body problem EBJ . In such a 



mapping the moment number corresponds to the number of particles interacting like charges 
on contact. Thus we speculate that the shoulders could be related to sudden changes in the 
character of the ground state as the particle number (moment number) increases. We will 
discuss this in more detail in the final section. 

For even smaller p values the curve starts to resemble the well known p = limit 
given by eq.(f|), and depicted as a flat line at In 2 in Fig. ^. From the figure one can 
graphically identify the value of (In J) as a function of p using the relation <iln( J^)/dn\ n ^Q = 
(In J). The quantity (In J) is a "free energy" that may reflect the sign transition. We have 
followed the value at intercept mentioned before as a function of p below p = 0.2. When the 
moments "zigzag" there are two possible extrapolations, while below the assumed transition 
the moments lead to a single prediction of the free energy. The results are depicted in Fig. 
[7|, where the curves merge around p c = 0.07 within the error of the extrapolation procedure. 
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Such a value coincides with our Monte Carlo prediction. 

One can validate the relevance of hierarchical lattices by checking the exact computation 



of the variable SJ/(J) with 5 J = J (J 2 ) — (J) 2 - Such a quantity has been discussed exten- 
sively in previous work Pf4|,|8|,f2~3}1 . As mentioned before, 5 J/ (J) was initially suggested as 
a candidate for a kind of order parameter that diverged exponentially above the transition 
and went to zero below. Observations by Shapir and Wang [[| showed, nevertheless, that 
path correlations (crossings) invalidated the vanishing of the parameter for any value of p. 
It has been argued that for small p there is a crossover from exponential growth (for p > p c ) 



to logarithmic growth (for p < p c ) flqj24}|. Shapir and Wang, on the other hand found a 
change from exp[|Zop(l — 2p)\2avL] for p < p c to exp[|Zog(2(l — 2p) 2 )\2L] for p > p c . Yet, 
they observe that the former result is incorrect because partial overlaps of pairs of walks 
should be accounted for. 

Simulations on regular lattices to date can only do very poorly in proving the surmised 
logarithmic behavior below p c . Here we have computed SJ/(J) to sizes L = 2 20 for various 
p values in a few CPU minutes. We have found a clear confirmation of logarithmic to 
exponential crossover as p increases. Figures [8] and [| show 5 J / (J) and its derivative as a 
function of L respectively. The scales used permit rapid identification of the corresponding 
behavior. It should be noted that, on euclidean lattices, it is reported the behavior reported 
is 5 J J (J) oc (lnL) M , where \i ~ 1 but depends weakly on p. 

On hierarchical lattices we can also demonstrate analytically that there is no transition 
in the variable 5J/(J). Following Cook and Derrida |HJ], eq.@, for the first two moments, 
can be written for general b as 

(Jm+i) = b(J m ) , 

(J 2 m+ i} = HJ 2 J 2 + Kb-l)(J m ) 4 . (6) 

Now, after defining J2(rn) = (J m ) 2 / '(</„) one can write a recursion relation for 5J / (J) as 



2 



oJm+1 \ I 



(J m+ i)y b 

n 



i-j1( m ) 



(7) 



It is simple to determine that j 2 has in general three fixed points, j 2 = 0, 1, 1/(6 — 1). For 
b > 2 (d e ff > 2) a critical fixed point arises and SJ/(J) exhibits a phase transition as NSS 
proposed. On the other hand, for b = 2 there are only two trivial fixed points; j 2 = 1 
is unstable and j 2 = is stable, indicating that SJ/(J) always diverges as found above. 
Values of j 2 close to one correspond to p — > 0, while j 2 close to zero correspond to p — > 1/2. 
Analyzing the behavior of the recursion for j 2 near the j 2 = fixed point one can derive 
from Eq.| that 5 J/ (J) ~ l/2exp(L(| In j 2 (0)\ + (l/2)j|(0))). The behavior close to j 2 = 1, 
which should be logarithmic, is also verified (numerically) although we have not arrived at 
a simple closed expression. In summary, hierarchical lattices provide similar results to those 
expected on euclidean lattices thus seeming a good testing ground for the sign transition. 

As a final word; we have computed higher order cummulants of J finding no features 
of special interest related to the transition. The only result worth mentioning is that 
hiiCj) 1 ^ / L = In 2 for p = 0, where Cj is the jth cummulant of J. In what follows we 
will take advantage of the special structure of hierarchical lattices to compute the full prob- 
ability distribution for J . 

V. PROBABILITY DISTRIBUTION FOR J 

Monte Carlo sampling of the distribution of J is handicapped by the models distribution 
broadness. For such reasons Wang et al || undertook an exact enumeration study to probe 
the NSS order parameter AP = P(J > 0) — P(J < 0). Because of the high computer 
demand of exact enumeration, they could only access sizes of L = 10 for all p. Here we use a 
scheme, on hierarchical lattices, permitting access to L = 16 exactly for all p and sampling 
of the distribution for L = 64. The procedure is as follows: As a hierarchical lattice is built 
recursively following a chosen motif, one can write the following recursion relation for the 
probability distribution. 

p w( J ) = II / Pi(jk)8(J - vim - V3V4)dvi, (8) 

i=i J -°° 
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where 771,2 and 773,4 denote contiguous elements on separate branches of the hierarchical 
lattice. Pi = p5(i] — 1) + (1 — p)S(i] + 1), where p is the sign probability discussed in previous 
sections. The number of possible outcomes for J or number of different paths goes as 2 2 ™ _1 
(32768 for L — 16, m = 5, and 2147483648 for L = 32, m = 6). This growth is extremely 
fast, although many J values will be degenerate for any particular disorder realization. Note 
that while L = 16 is easily accessible, going an order further, puts the calculation out of 
reach, no intermediate sizes being available on hierarchical lattices. For L = 32 we have 
resorted to a coarse-graining procedure in the following manner: the exact results for L = 16 
involve 175 terms which we cannot evolve exactly to the next order. Nevertheless, we can 
make a coarse-grained distribution by averaging J occurrences in groups of 7 to obtain 25 
different values. One can then go to up to L = 64 by repeating this procedure. Beyond such 
a size, the coarse graining procedure does not incorporate sufficient detail to see anymore 
changes in the distribution, so within our resolution we have achieved its limit form. 

Fig. [I~C)| a shows the probability distribution for L = 16 for significant values of p. The 
probability distribution is astonishingly complex, even for small sizes, revealing rich inter- 
ferences in the paths sums J. Note that it would be hopeless to sample the distribution 
P(J) using Monte Carlo as there are sixty to one hundred and thirty orders of magnitude of 
probability. Fig. |i~0l b shows different p values for a sample L = 32 as a function of p using 
the coarse-graining procedure described above. As expected the distribution is symmetric 
for p = 0.5 and gains asymmetry (AP 7^ 0) as p moves towards zero. Note that P(J) falls 
slower than exponential, on average, about the peak value. It is notable the speed with 
which the distribution appears symmetric beyond p = 0.1. This feature is understood in the 
'zigzag' behavior of the moments, where odd moments are much smaller than even moments 
and their separation increases exponentially as p is increased. 

Having the information of the exact distribution one is also able to obtain the exact order 
parameter AP introduced in section III. No qualitative differences were found with curves 
reported in figures 1 and 2 at least to sizes L = 32, so sampling of AP, involved in Monte 
Carlo, seems to be good enough to draw the conclusions about the transition (see section 
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III). 

In Fig. |TT] we have depicted the distribution for L = 64 and p = 0.5 without joining 
the points for the probability amplitude (as was done in Fig. |H]). A fractal structure is 
apparent; The whole distribution, in the shape of an approximate triangle, is built from, 
scaled, identical triangular structures up to the resolution achieved by the coarse-graining 
procedure. Similar complexity is expected for the sign problem on Euclidean lattices. 

An interesting final point to make in this section is that, in view of the unique relation 
between distribution and moments (see section IV), it is possible to use known inversion 
formulas [f^] . In this way one could derive the limiting distribution exactly to any order 
desired. 

VI. SUMMARY AND DISCUSSION 

We have provided evidence for the existence of a phase transition for the directed path 
sign model on hierarchical lattices. Nontrivial finite size scaling of the order parameter close 
to the transition, points to a second order phase transition as found from mean field type 
approaches. From numerical computations, the threshold on diamond hierarchical lattices 
is p c = 0.071 ± 0.001 and the correlation length exponent is v = 1.85 ± 0.07. The latter 
exponent is very different from that of percolation on the same lattice v = 1.635... 

The study of exact moment recursion relations for (J n ) led us to the definitive conclusion 
that the moments uniquely determine the probability distribution, according to Carleman's 
theorem [p2 ]. Using extrapolations of the derivative of integer moments (d(J) n /dn) to 



n = 0, we were able to find a "free energy" (In J). Such a free energy splits into two 
possible extrapolations (from even and odd moments) as one goes through the transition 
point by increasing p. The latter transition point coincides with that found in Monte Carlo 
simulations of the sign transition. We have not completely interpreted this connection in 
the present paper. Furthermore, we studied the high moments of the partition function J 
below the transition, and found a very interesting non-monotonic behavior including step 
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structures that propagate on the moment number axis, as p changes. 

Using the fact that moments can be computed exactly we studied the celebrated ratio 
SJ/(J) proposed by NSS. We have shown, analytically that indeed in d e ff = 2 the ratio does 
not show a transition as suspected numerically P,|8|,|2"4|| on regular lattices. Furthermore, we 
have shown that hierarchical lattices exhibit the same logarithmic to exponential crossover 
for SJ/{J) surmised in references || and [p4 ]. 



Finally, we studied a recursion relation for the full probability distribution for J, finding 
an extremely complex structure even for systems as small as L = 16. Previous remarks by 
Roux and Coniglio |L8| of anomalous accumulation of probability at J = are confirmed. 
Nevertheless, their claim that the hierarchical lattice becomes essentially one dimensional 
for large L, and thus, the probability distribution should approach a Gaussian is not borne 
out from our results. One obvious difference is that for a Gaussian all cumulants, larger than 
two, should be zero which is in disagreement with exact results of section IV. No evident 
signal of the transition, beyond that already inferred from the order parameter AP, is found 
from the full probability distribution. 

Medina and Kardar [[3] have studied the moments for the sign problem interpreting them 
as partition functions for n-body one dimensional Hamiltonians with contact interactions. 
Most of the focus, however, has been on the low n behavior that yields cummulants of In J . 
Nevertheless, it would be interesting to interpret the findings of this paper, regarding high 
moments, in the light of a many body-theory A previous effort by Zhang |23| focused on the 



Hartree-Fock approximation valid only for a large number of particles (higher moments). In 
Zhang's approach the sign model was equivalent to finding the ground state of the many 
body Hamiltonian 

2n \ 

E d i + E wA^i - Xj) *(xi...x„) = E (n)*(x 1 ...x n ) (9) 

i=l i>j ) 

where e, is a charge that acts via contact interaction of the ith particle: ei = 1 for 1 < i < n 
and Cj = —1 for n < i < 2n. Zhang's approach yields E$ oc n 2 . Our findings predict, 
from the relation ln(J n )/L = E , E = ^n for large n (see also ||), where 7 increases as 
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p — > 0. For lower n the behavior is non-trivial and is certainly not represented a simple 
powerlaw of n. Therefore, Zhangs results represent some kind of intermediate regime. A 
more detailed solution of eq. (|9|) might yield the 'solitonic' patterns reported here (see 
section IV) which are not well understood. As speculated in section IV, the ground state 
formed by particles with attractive and repulsive interactions might change, suddenly, at 
critical particle numbers generating discontinuities in the derivative of ln(J n ). More work 
is needed in this direction. 

The highly non-monotonic behavior displayed by the moments calls for caution regarding 
the regime of validity of moments dependencies on the moment number n reported in the 
literature |]3|,|23[|. Claims of a non-unique relation between moments and the probability 
distribution were based on expressions only valid in the n — > limit, which is clearly 
unrelated to the constraints of Carlemans theorem [25|| . Obviously, the conclusions of this 



paper are only valid in the measure to which hierarchical lattices agree with continuum 



results. For a discussion of the latter point see reference [13 
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FIGURES 
FIG. 1. Hierarchical lattices are built by repeating a chosen motif; each bond turns into a 

diamond recursively. The figure shows successive iterations of the lattice and the corresponding 

length L between end points i and /. Examples of a directed path at each order are indicated by 

contiguous arrows. 

FIG. 2. The figure shows the order parameter AP = P{J > 0) — P{J < 0) as a function of 
the sign probability p for system sizes indicated. Averages were performed over more than 20000 
realizations of randomness. Note the formation of a plateau at AP = 1 for small p. 

FIG. 3. Same data as in Fig. ^ after collapsing the curves for different system sizes. The 
appropriate choices for p c and u, the transition threshold and the correlation length exponent, are 
indicated. 

FIG. 4. The figure shows the value of p c (L), evaluated from the peaks of the derivative of the 
order parameter, as a function of (l/L) 1 '". The last five sizes from L = 32 to L = 512 have been 
fitted by least squares to yield the asymptotic value p c (oo) = 0.072 indicated. 

FIG. 5. The moments ln{J n )/L as a function of the moment number n for the lattice sizes 
indicated. The figure shows the rapid convergence to an asymptotic result. In the inset, we show 
that while the initial moments grow faster than exponential they nevertheless grow slower than n\, 
so there is a unique relation between moments and probability distribution. 

FIG. 6. The derivative of ln(J n )/L as a function of the moment number, for size L = 2 15 , (last 
six orders collapse onto the same curves.) As p decreases the curves approach the asymptotic value 
In 2. Note the change from the sawtooth behavior above p = 0.07 to collinearity. Shoulder features, 
developing at higher n, move almost undeformed in the positive n direction as p increases. 

FIG. 7. The "free energy" (In J)/L, as discussed in the text, as a function of p. Note that the 
curves extrapolated from even and odd moments merge around the threshold for the transition 
obtained from Monte Carlo. 
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FIG. 8. The figure shows the behavior of 8 J/ (J) as a function of L for the p values indicated. 
As the plot is semi-log the exponential behavior above p = 0.03 is evident. 

FIG. 9. The figure shows the derivative of the data in Fig. @. Here the logarithmic behavior of 
8J / (J) is evident. The dotted line is a guide for the eye for 1/L behavior. 

FIG. 10. a) The exact probability distribution for L = 16 at the p values indicated. Note the 
self-affine structure; the central peak is repeated at subsidiary local maxima, b) The figure shows 
the coarse-grained distribution for L = 32. 

FIG. 11. The probability distribution for L = 64 and p = 0.5 using the coarse-grained distri- 
bution for L = 32. To emphasize the self-similar structure we have not to joined the points in the 
graph as in Fig. [ji]. 
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